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ABSTRACT 



We present a fitting formula for the luminosity and angular diameter 
distances in cosmological models with pressureless matter, a cosmological 
constant and zero spatial curvature. The formula has a relative error of less 
than 0.4% for 0.2 < f2 < 1 for any redshift, and a global relative error of less 
than 4% for any choice of parameters. 



Subject headings: Cosmology: theory - distance scale 
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1. Introduction 

It is presently fashionable to consider spatially flat cosmological models with a 
cosmological constant (Ostriker and Steinhardt 1995). Several projects are underway to 
measure the magnitude of the cosmological constant, usually based on the luminosity- 
distance (Perlmutter et al. 1997), the angular-diameter distance (Pen 1997), or 
volume-distance (Kochanek 1996). Unfortunately, these distances are only expressible in 
terms of Elliptic functions (Eisenstein 1997). In order to simplify the repeated computation 
of difficult transcendental functions or numerical integrals, we present a fitting formula with 
the following properties: 

1. it is exact for f2 — ► 1~ and fi — ► + at all redshifts. 

2. The relative error tends to zero as z — > oo for any value of Q . 

3. For the range 0.2 < O < 1, the relative error is less than 0.4%. 

4. For any choice of parameters, the relative error is always less than 4%. 

Without further ado, the luminosity distance is given as: 
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r)(a,n ) = 2Vs 3 + 1 



* 3 ^ ^ w 

We have used the Hubble constant H , and the pressureless matter content f2 . We recall 
that an object of luminosity L has flux F = L/(And L ). 
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2. Approximation 

Spatial flatness allows us to rewrite the Friedman-Robertson- Walker metric as a 
conformally flat spacetime 

ds 2 = a 2 {-drf + dr 2 + r 2 dQ). (2) 

7] is the conformal time, and r is the comoving distance. We have set the speed of light 
c = 1 . The scale factor a can be normalized in terms of the redshift a = 1/(1 + z) . The 
Friedman equations determine 



' da\ 2 



aQ + a 4 {l\ (3) 



\drj J 

where spatial flatness requires Qq + Q\ = 1. A change of variables to u = as where 
s 3 = (1 — Qq)/Qq allows us to express (§) parameter free 



's 3 + l f as du 

(4) 



u 



We can asymptotically approximate (f|) as 



\u- n + n(2X) 2n+1 u- 1 + (2X) 2n l 2 " (5) 



where X = (J °° du/y/u 4 + u) _1 = 3v / 7r/(r[|]r[|]) ~ 0.3566. n is a free parameter, and at 
this stage a choice of n = 3 approximates all distances to better than 14% relative accuracy. 
Equation (|5|) satisfies the following conditions: 1. it converges to (H) as u — > and u — >• oo. 
2. Its derivative converges to the derivative of rj as u — > oo. To improve on (|5|) we consider a 
polynomial expansion of u~ 1 in the denominator with two free parameters by setting n = 4: 



y' f!±l + ClU -3 + C2U -2 + 4( 2X) 9 n- 1 + (2X) 8 ] _1/8 . (6) 

We will now choose the coefficients Ci, C2 to minimize the relative error in the approximate 
luminosity distance di, 

eoo = WiL-d^/d^U (7) 
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where the subscript oo indicates the infinity norm, i.e. the maximal value over the domain. 
The error in (0) tends to be dominated by z = (see for example Figure |]), for which we 
can express 

\At 4 + u-r- - 1 
du 

Globally optimizing (H) allows us to reduce the error to about 2%. But we choose the 
following trade-off for current cosmological parameters: We want to minimize (0) over the 
range 0.2 < Qq < 1, which covers the popularly considered parameter space. In that range, 
we find through a non-linear equation solver that c\ = —0.1540 and C2 = 0.4304. This 
allows us to trade the global error of 2% for a global error of 4%, while reducing the error 
in the range of interest to 0.4%. The global error surface plot is shown in Figure The 
error at z —>■ is shown in Figure |I[ At small z, any errors in fj are amplified, even though 
the global errors in r)(z) are generally significantly smaller. Figure (0) shows the fit for the 
conformal time fj and its residual, which is accurate to 0.2% globally. 

One can always express the luminosity distance as a power series expansion around 
z = using Equations ([[]) and (U). For a flat universe, one obtains 

-' , f i 3^ \ 2 , irA^° 3 



-fd L = z+\l- -fioj z 2 + (9fi - 10)y z A + ■■■ (9) 

The series converges only slowly for z ~ 1: even when expanded to sixth order in z, the 
maximal relative error in the interval < z < 1 and 0.2 < Qq < 1 using the series expansion 
in (|) is 37%. 



3. Conclusion 



We have presented a simple algebraic approximation to the luminosity distance di and 
the proper angular diameter distance dA = d^jiX + z) 2 in a flat universe with pressureless 
matter and a cosmological constant. 
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z -> 




Fig. 1. — Relative error at z — > 0, which dominates the global error budget. The minmax 
error was optimized in the range 0.2 < Q < 1, which is the area to the right of the vertical 
dashed curve. The error goes to zero as a powerlaw at both ends of the graph, which is 
shown in the two insets. 
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Fig. 2. — Global error surface. Our minmax fit minimized all errors for log 10 > —-7, 
where the error is less than 0.4% for all redshifts. 
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Fig. 3. — Exact and approximate conformal time. The lower figure shows the conformal 
time as a function of u — [(1 — VL$) /VLq] 1 ^ / {1 + z) for the exact integral (solid line) and the 
approximation of equation ([T]) (dashed). These two lines are practically indiscernible. The 
upper figure shows the relative error (fj — rj)/r}. 



